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Abstract. The number of constituent quarks (NCQ-) scaling of hadrons and the number of constituent 
nucleons (NCN-) scaling of light nuclei are proposed for nuclear modification factors ( R cp ) of hadrons and 
light nuclei, respectively, according to the experimental investigations in relativistic heavy-ion collisions. 
Based on coalescence mechanism the scalings are performed for pions and protons in quark level, and 
light nuclei d(d) and 3 He for nucleonic level, respectively, formed in Au + Au and Pb + Pb collisions and 
nice scaling behaviour emerges. NCQ or NCN scaling law of R cp can be respectively taken as a probe for 
quark or nucleon coalescence mechanism for the formation of hadron or light nuclei in relativistic heavy-ion 
collisions. 

PACS. PACS-key 25.75.Gz, 12.38.Mh, 24.85.+p 

1 Introduction 

Quark-gluon plasma (QGP) which was predicted by the 
Quantum chromodynamics (QCD) [T] can be formed in 
ultra-relativistic heavy-ion collisions, at the Relativistic 
Heavy-Ion Collider (RHIC), Brookhaven National Labo¬ 
ratory [5JEM1II3 and at the Large Hadron Collider (LHC), 

CERN, and the energy dependence of quark matter prop¬ 
erties is thought as one of important observables to study 
various aspects of the QCD phase diagram mmwmm 
n~21fITTirT~TlH771 in the beam energy scan program at RHIC 
|16j . Extensive experimental results in Au + Au collisions, 
including particle and anti-particle production and spec¬ 
tra, the nuclear modification factors (NMF) Raa or R cp 
of hadrons and baryon-to-meson ratios at y /snn = 200 
GeV/c have been reported p[51l¥l[5I IT71IT51[roiBIJl[iT] . Some 
thermal parameters of the QGP are extracted, for instance 
the chemical freeze-out temperature T c h , the baryon chem¬ 
ical potential pb and the strangeness chemical potential 
ps etc mm, which helps us to understand the phase 
structure at high temperature. 

In the viewpoint of experimental probes, the nuclear 
modification factor (NMF) R cp is a very useful observ¬ 
able for studying the feature of QGP (331I3H1 . Extensive 
studies have been carried out both experimentally and 
theoretically. These studies indicate that the NMF, which 
can be represented either by the modification factor be- 
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tween nucleus-nucleus (AA) collisions and proton-proton 
(pp) collision Raa or the one between the central col¬ 
lisions and peripheral collisions R cp , can give quantita¬ 
tive properties of the nuclear medium response when the 
high speed jet transverses it. In high transverse momen¬ 
tum (pp ) region, NMF is suppressed owing to jet quench¬ 
ing effect in hot-dense matter and thus has become one of 
the robust evidences on the existence of the Quark-Gluon- 

dial flow boosts or the Cronin effect GEEM31 competes with 
the quenching effect and enhances the NMF, which has 
been also demonstrated by the Relativistic Heavy-Ion Col¬ 
lider (RHIC) beam energy scan (BES) project 33]. Re¬ 
cently, the study of nuclear modification factor was for 
the first time extended to intermediate-energy heavy-ion 
collisions and the radial flow effect on proton NMF has 
been quantitatively investigated by Ma’s group [841185] . 


It will be of very interesting to check if the R cp of n 
and p can be scaled together by a number of constituent 
quark (NCQ) scaling law as elliptic flow’s NCQ-scaling [B] 
mm- On the other hand, R cp of light nuclei is not yet 
addressed in relativistic heavy ion collisions. In the present 
study, apparent differences between the R cp of light nuclei 
are found for different nucleon number. With the help of 
the coalescence model [37], we also found a way to fit them 
by a number of constituent nucleon (NCN) scaling. In the 
following texts, we will discuss these phenomena. 
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2 Number of constituent quark (NCQ) 
scaling of hadron’s R cp 


In relativistic heavy ion collisions the production of hadrons 
at intermediate transverse momentum can be described by 
coalescence mechanism '28], i.e. 


p d 3 N h 

h d*P h ~ Bn 


P d 3 N q 
q d?P a 


(1) 


where n q is the number of constituent quarks of a hadron 
and the coefficient B nq is the probabilities for n q quarks 
to hadron coalescence. From the coalescence mechanism 
it was found that elliptic flow of hadrons can be scaled 
to the number of constituent quark msm- And the 
nuclear modification factor R cp is defined as mm, 

p , [d 2 N/p T dydp T /(N bin )} cent ™ 1 

cp[PT> [d?N/p T dydpTl{N bin )]P eri P heraV 1 



P t (GeV/c) 


Fig. 1. (Color online) R cp for pions and protons in Au+Au 
collision at y/ snn = 200 GeV. Original data of R cp are taken 
from ref. m- 


where ( Nb % n ) is the average number of binary nucleon- 
nucleon collisions per event. The ratio of (Nbi n ) central to 
(Nbin) perzpharal can be found in data .17] for Au+Au colli¬ 
sions at y/SNN = 200 GeV and in data [+711+01 for Pb+Pb 
collisions at y/S nn = 2760 GeV. The jets and high pt 
particles created in the early stage will lose most energy 
through interactions in the evolution of the system in cen¬ 
tral A+A collisions and it is confirmed by suppression of 
R cp at high pr in experiments in relativistic heavy-ion 
collisions. Figure |T] displays R cp for pions and protons in 
three central-peripheral pairs, where the original data are 
from ref. |41j . Of course, obvious differences between pi¬ 
ons and protons are observed. In this paper, however, we 
will focus on another topic, i.e. if there is number of con¬ 
stituent quark scaling of R cp for hadrons as that for ellip¬ 
tic flow at intermediate transverse momentum. 

From eq. HD and eq. m, it can be deduced, 


KM) 


R 


n q , central 


B 


n q ,peripheral 


l/n q 

(R cp (n q -pr)) 1 ng 


( {N bi n) c \ 1/nq 1 

\(N bin )p) 


( 3 ) 


R* is the scaled nuclear modification factor to scale me- 

'-'P 

son’s and baryon’s R cp together. Unfortunately the coef¬ 
ficient B nq can not be determined experimentally. So we 
can try to scale hadron’s R cp by the following formula, 


KM) = ( Rcp{n q -p T )) 


1/riq ( (NbinY 


\{Nbi„)P 


l/n q -l 


KM) = 


K p (pt) 


B n 


i q peripheral 
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( 4 ) 


Figure [2] presents R* p of pion and proton as a function 
of transverse momentum in Au+Au collisions at y/spjN = 
200 GeV, where the original data of R cp are from ref. SB- 
The value of R* p of proton is scaled by a factor listed in the 
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Fig. 2. (Color online) Number of constituent quark scaling of 
Rcp for pion and proton in Au+Au collision at y/ snn = 200 
GeV. Original data of R cp are taken from ref. mi- 


fig - ei which can scale pion and proton R* p together. It can 
be seen from the listed factor which shows the difference 
between pion’s R* p and protons is a constant factor with 
centrality dependence. 

Similarly, R cp of pions and protons for the Pb+Pb col¬ 
lisions at y/SNN = 2760 GeV are plotted as a function of 
Pt in fig. [3] for different central-peripheral pair combina¬ 
tion. Again, significant differences are observed between 
pions and protons. After the number of constituent quark 
scaling for R cpi we plot fig. Q] for demonstration of R* p 
of pion and proton, where the original data of R cp are 
from [33Mi]E2M3] . The scaled factor of R* p for proton is 
increasing with the increasing of the centrality of numer¬ 
ator of R* p , which takes similar centrality dependence in 
LHC as in RHIC energy. So it implies that the number of 
constituent quark scaling of R cp for hadron supports the 
viewpoint of quark coalescence mechanism for hadron for¬ 
mation in the intermediate transverse momentum range 
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Fig. 3. (Color online) R cp for pion and proton in Pb+Pb collision at y/ snn = 2760 GeV. Original data of R cp are taken from 
refs. fggU40U42U43) . 
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Fig. 4. (Color online) Number of constituent quark scaling of R cp for pion and proton in Pb+Pb collision at y/ snn = 2760 
GeV. Original data of R cp are taken from refs. I39II40II42II43I . 


for relativistic heavy-ion collisions and the difference of 
the scaled R cp between pions and protons keeps a con¬ 
stant factor which displays a centrality dependence. 

To summarise this section, we find that the constituent 
quark scaling holds for the nuclear modification factor of 
hadrons. It indicates the energy loss of hadrons is essen¬ 
tially originated from the partonic stage which leads to 
the same energy loss factor per constituent quark, which 
is consistent with the concept of quark-gluon plasma. 


3 Number of constituent nucleon (NCN) 
scaling for R cp of light nuclei 


While the collision system reaches a kinetic freeze-out 
stage, nucleons could be coalesced to light nuclei. In this 
case, the coalescence mechanism for the formation of light 
nuclei is similar to formula m with considering the same 
distribution of protons and neutrons, 


F d 3 N A 

Ea 1hr- Ba 


p d 3 N P 
P d?P p 


E, 


d 3 Nr, 


A-Z 


=+*++)"■ < 5 > 

where 

P P = Pn = Pa/'A, (6) 

and the coefficient B A ~ (1/V")^” 1 is related to the fire¬ 
ball volume in coordinate space and depends on momen¬ 
tum mmm- The coefficient B A can be extracted from 
data or calculated by coalescence mechanism HUBS] . The 
number of constituent nucleon (NCN) scaling of light nu¬ 
clei’s R cp can be obtained through replacing n q by A in 
eq. (£[]>. 

If there exists a similar scaling of R cp for light nu¬ 
clei as that for elliptic flow [35M71I38] . it certainly sup¬ 
ports the coalescence mechanism for formation of light 
nuclei at kinetic freeze-out stage. Our earlier prediction 
on number-of-nucleon scaling of elliptic flow mm was 
recently supported by the STAR data from Au + Au colli¬ 
sions [4$J . Even though the energy domain is huge different 
between our prediction mm and the STAR data gSJ, 
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Fig. 5. (Color online) R cp of light nuclei of (a) Pb + Pb at 17.2 GeV, (b) Au + Au at 200 GeV, and (c) Pb + Pb at 2760 
GeV, respectively . Original data are taken from the NA49 Collaboration m for Pb+Pb collisions at y /snn = 17.2 GeV, the 
PHENIX Collaboration I52II53I for Au+Au collisions at y/SNN = 200 GeV and the ALICE collaboration I.39ll40ll421l43ll54l for 
Pb+Pb collisions at y/SNN = 2760 GeV. 
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Fig. 6. (Color online) Number of constituent nucleon (NCN) scalings for R cp of light nuclei of (a) Pb + Pb at 17.2 GeV, (b) 
Au + Au at 200 GeV, and (c) Pb + Pb at 2760 GeV, respectively . Original data are taken from the NA49 Collaboration [SI] 
for Pb+Pb collisions at yj snn = 17.2 GeV, the PHENIX Collaboration 152 11531 for Au+Au collisions at y/SNN = 200 GeV and 
the ALICE collaboration I39II401I4~2U43U54| for Pb+Pb collisions at y/SNN = 2760 GeV. 


the scaling phenomenon is the same, which indicates that 
the production of light nuclei stems from the nucleonic co¬ 
alescence mechanism rather than quark coalescence even 
in ultra-relativistic energy. The same concept of nurnber- 
of-nucleon scaling of elliptic flow was also followed by a 
dynamical model |3S] as well as a multiphase transport 
(AMPT) model [501 . 

Here we consider the coalescence mechanism for light 
nuclei by Eq. ©■ As the yields of p(p) and n(n ) are simi¬ 
lar, we have the yields of light nuclei with A-mass number 
obey the Eq. ©■ Before the scaling, we show the R cp of 
light nuclei of Pb + Pb at 17.2 GeV [51], Au + Au at 200 
GeV [521155] , and Pb + Pb at 2760 GeV [551I521H51I551 
HD] , respectively. Again, the curves for different light nuclei 
are not collapsed together. After the scaling for nucleon 
number, fig. [6] shows the number of constituent nucleon 
on (NCN) scaling for R cp of light nuclei, R* p . Figure[6](a) 

represents the R* p of proton and deuteron in Pb + Pb col¬ 
lisions at y/SNN = 17.2 GeV and the data are taken from 


the NA49 collaboration [51]. Figure [6] (b) gives the R* p of 
(anti-)proton and (anti)deuteron in Au+Au collisions at 
y/SNN = 200 GeV and the data are from PHENIX col¬ 
laboration [52U53j . And fig. [6] (c) shows the R* p of proton, 
deuteron and 3 He in Pb + Pb collisions at tJsnn = 2760 
GeV and the data are from ALICE collaboration [3911401 
mmm- From these plots, it was found the nuclear mod¬ 
ification factor (R cp ) of deuteron and 3 He can be scaled to 
proton’s R cp after performing the number of constituent 
nucleon (NCN) scaling. There is a constant factor for R* p 
of light nuclei to that of protons. Although the NCN scal¬ 
ing for R cp of light nuclei is studied at three different col¬ 
lision energy y/SNN from 17.2 GeV to 2760 GeV, there 
needs more experimental results of energy and centrality 
dependence of light nuclei production to explore the con¬ 
stant factor for R* p of light nuclei which is related to the 
coalescence parameter from eq. ©■ 

In light of this study, we can investigate the R* p for 
light nuclei experimentally to distinguish their formation 
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mechanism in heavy ion collisions. Current data support 
that the (anti-) light nuclei are coalesced by the (anti-) 
nucleons. In fact, the constituent nucleon number scaling 
of elliptic flow for light nuclei which was firstly proposed 
in intermediate-energy heavy-ion collisions 
also supports the coalescence mechanism for the formation 
of light nuclei at kinetic freeze-out stage in the reaction 
system. Actually, at RHIC, the elliptic flow of (anti-) light 
nuclei data of Au + Au at 200 GeV [38] gives the evidence 
of the nucleon-number scaling as proposed by Ma’s group 
[38] . Furthermore, the observation of anti-hypertriton m 
and anti-helium 4 m at STAR also indicates that the 
anti-nucleon coalescence is responsible for their produc¬ 
tion yield in the basis of attractive interaction between 
anti-nucleons pTlfHR] . 

To summarise this section, we find that the constituent 
quark scaling holds for the nuclear modification factor of 
light nuclei in relativistic heavy-ion collision. It indicates 
the energy modification of light nuclei is essentially origi¬ 
nated from the hadronic stage rather than protonic stage, 
which leads to the same energy loss factor per constituent 
nucleons. On the other hand, in contrast to the scaling be¬ 
havior of R cp for hadrons, the scaled R cp displays a rising 
dependence as a function of pr, which could stem from 
the larger radial flow or multiple nucleon scattering for 
(anti-) light nuclei as demonstrated in refs. [3311851 . 

4 Conclusion 

We propose, for the first time, that the number of con¬ 
stituent quark (NCQ-) scaling of hadrons and the number 
of constituent nucleons (NCN-) scaling of light nuclei for 
nuclear modification factor R cp , respectively, in this work. 
The NCQ-scaling of R cp of hadrons was applied to the re¬ 
sults in Au + Au collisions at yjS nn = 200 GeV from the 
RHIC-STAR and in Pb+Pb collisions at yj snn = 2760 
GeV from the LHC-ALICE. The scaled R cp presents a 
constant difference between pions and protons at interme¬ 
diate transverse momentum and the constant difference 
factor shows a centrality dependence trend. This is a sup¬ 
plement for number of constituent quark scaling besides 
elliptic flow which supports the coalescence mechanism for 
formation of hadrons. It is also consistent with the forma¬ 
tion of quark-gluon plasma in Au+Au collisions at y/SNN 
= 200 GeV. 

Based on the coalescence model, the scaling of R cp has 
also been tested for light nuclei where the number of con¬ 
stituent nucleon (NCN-) scaling is taken into account. The 
results display that the R cp spectra of different light nu¬ 
clei become the one as protons after the scaling. This im¬ 
plies that the scaling factor for R cp of light nuclei should 
be a constant value and would not be distorted for the 
structure of R cp . This scaling behaviour can be taken as a 
nucleon coalescence mechanism of light nuclei, and illus¬ 
trates that those (anti-) light nuclei are originated from 
the hadronic stage. 
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